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In this Article we review some recent progresses in the field of non-equilibrium linear response the- 
ory. We show how a generalization of the fluctuation-dissipation theorem can be derived for Markov 
processes, and discuss the Cugliandolo-Kurchan '3] fluctuation dissipation relation for aging systems 
and the theorem by Franz et. al. ^] relating static and dynamic properties. We than specialize the 
subject to phase-ordering systems examining the scaling properties of the linear response function 
and how these are determined by the behavior of topological defects. We discuss how the connection 
between statics and dynamics can be violated in these systems at the lower critical dimension or as 
due to stochastic instability. 
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I. INTRODUCTION 

The fluctuation-dissipation theorem [l[ (FDT) is one of the fundamental accomphshments of hnear response theory 
appUed to equilibrium systems. According to the FDT a response function x, describing the effects of a small 
perturbation exerted on a system, is linearly related, via the equilibrium temperature T, to a correlation function C 
of the the system in the absence of the perturbation. In the language of magnetic systems, which we shall adopt in the 
following, one usually considers the application of an external magnetic field /i, and x is the magnetic susceptibility. 

In recent years there has been a considerable interest, arisen in different fields such as turbulent fluids 0, disordered, 
glassy [3, '3,^ and aging systems in the generalization of the results of linear response theory to out of equilibrium 
systems. Differently from equilibrium statistical mechanics, where a well funded and controlled theory is available, 
there is not nowadays a theorem of a generality comparable to the FDT for non-equilibrium states. Nevertheless, some 
interesting progresses have been done in understanding some particular aspects of non-equilibrium linear response 
theory, some of which will be discussed in this paper. 

A flrst basic question regards the possibility of generalizing the FDT. Namely, the question is whether also away 
from equilibrium it is still possible to relate the response function to properties of the unperturbed dynamics, possibly 
in the form of correlation functions. A positive answer to this question exists when the time evolution is Markovian 
and described by a differential equation of the Langevin type Q or for systems described by a master equation . 
In this case, the response function is related to correlation functions of the unperturbed systems, which, however, 
are not only the correlation C involved in equilibrium. These results are particularly important since they allow the 
study of the response function without considering the perturbed system, which is generally more complicated. 

Once a relation between response and correlations is established, at least in the restricted framework of Markov 
processes, the natural question is which piece of information, if any, can be learned from it about the non-equilibrium 
state. In equilibrium the linear relation between x and C is universal, and the coefficient dx/dC entering this relation 
is T. In the restricted area of aging systems it has been shown that the x(C') relation still bears an universal character, 
although weaker than in equilibrium. This is because the theorem by Franz, Mezard, Parisi and Peliti connects 
x(C) to the equilibrium probability distribution of the overlaps Peqio) and different statistical mechanical systems 
can be classified into few universality classes on the basis of their P^q (q) according to the replica symmetry breaking 
character of the ground state [13 . Moreover, dx/dC can be interpreted as an effective temperature [J]. 

These results promoted linear response theory as an important tool to investigate the non-equilibrium behavior or 
even the structure of equilibrium states of complex systems, such as spin glasses, which are hard to equilibrate, where 
Peqiq) can be better inferred from a non-equilibrium measurement of xi^)- However, in order for these studies to be 
sound, a basic understanding of the out of equilibrium behavior of the response function is required. Instead, already 
at the level of coarsening systems, which can be considered as the simplest paradigm of aging phenomenon, where a 
satisfactory general analytic description can nowadays be given by means of exactly solvable models or approximate 
theories, the scaling properties of the response function are non trivial and still far from being understood. Notably, 
the statics-dynamics connection stated by the theorem [9] is not always fulfilled in coarsening systems. 

In this paper we review some recent progresses in the field of non-equilibrium linear response theory. The focus is 
mainly on aging systems and, in particular, on phase-ordering, which, because of its relative simplicity, is better suited 
for a thorough analysis. The article is organized as follows: Sees. [Ill IIIIl IIVI and |V] are of a general character; here 
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we fix up the basic definitions, discuss the generahzation of the FDT for Markov processes, introduce the fiuctuation- 
dissipation relation and review the theorem Q which hnks statics to dynamics. In Sec. IVII some of the concepts 
introduced insofar are speciahzed to the case of phase-ordering kinetics. After a general description of the dynamics 
in Sec. IVI Al the behavior of the response function is reviewed in Sec. IVIBI In particular, in Sec. IVIB II the scaling 
properties of x are discussed and in Sec. IVIB 21 it is shown how, in the case of a scalar order parameter, the exponents 
can be related to the roughening properties of the interfaces. Sec. IVI CI contains a discussion of how the connection 
between statics and dynamics is realized or violated in coarsening systems. Some open problems are enumerated in 
Sec. I VIII and the conclusions are drawn. 



II. BASIC DEFINITIONS 



Let us consider a system described by the Hamiltonian Hq. The autocorrelation function of a generic observable O 
between the two times s and t > s is 

cit,.s)^m)ois)), (1) 

where (. . .) is an ensemble average. Switching on an impulsive perturbation h(s) at time s which changes the 
Hamiltonian Ti ^ Ti -\- ATi — H — hO, the linear (impulsive) response function is given by 



dh{s) 

The integrated response function, or dynamic susceptibility, is 



(2) 

h=0 



t 



Xit,s)= / R{t,t')dt' (3) 



and corresponds to the response to a perturbation switched on from s onwards. 

In equilibrium, time translation invariance (TTI) holds, so that all the two time quantities introduced above depend 
only on the time difference t = t — s. The FDT reads 

TRir)^-'-^, (4) 

where T is the temperature, or, equivalently, for the integrated response 

Tx{t) - C{r = 0) - C{r). (5) 



III. OFF-EQUILIBRIUM GENERALIZATION OF THE FLUCTUATION DISSIPATION THEOREM 

FOR MARKOV PROCESSES 

Consider a system with an order parameter field (f>{x) evolving with the Langevin equation of motion 

d(t){x, t) 



dt 



B[q^{x,t)]+r^{x,t) (6) 



where B [(j){x, t)] is the deterministic force and r](x, t) is a white, zero-mean Gaussian noise. In this framework a 
generalization of the FDT was derived in Let us recall the basic elements, referring to [3] for further details. From 
Eq. ([S]), the linear response function is simply computed as the correlation function of the order parameter with the 
noise 

2TRit,s) = {^{x,t)rjix,s)) (7) 

where T is the temperature of the thermal bath and i > s by causality. It is straightforward to recast the above 
relation in the form 
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where 

C{t, s) = t)<P{x, s)) - , t)) , s)) (9) 

and 

A{t, s) = i {(<^(f , t)B , s)]) - {B [<^(f , t)] 0(f , s))} (10) 

is the so called asymmetry. Eq. (jS]), or ([7]), qualifies as an extension of the FDT out of equilibrium, since in the right 
hand side there appear unperturbed correlation functions. When time translation and time inversion invariance hold, 
so that A{t, s) = and dC{t, s)/dt = —dC{t, s)/ds, it reduces to the equilibrium FDT (j4]). Let us mention that this 
equation holds [8| in the same form both for conserved order parameter (COP) and non conserved order parameter 
(NCOP) dynamics 

The next interesting question is whether one can do the same also in the case of discrete spin variables, where the 
kinetics is described by a master-equation, there is no stochastic differential equation and, therefore, Eq. ((7]) is not 
available. A first approach to this problem was undertaken in Refs. [13, [13, EB, EB] where a relation between the 
response function and particular correlators was obtained. As we shall discuss briefly below, however, their results 
cannot be qualified as generalizations of the fluctuation-dissipation theorem. Instead, in what follows we scketch how 
(details can be found in ,8,]), an off- equilibrium generalization of the FDT, which takes exactly the same form as 
Eqs. (|8ll0p and which holds, as in the Langcvin case, for NCOP (spin flip) and COP (spin exchange) dynamics can 
be derived also in this case. 

Let us consider a system of Ising spins ai — ±1 executing a Markovian stochastic process. The generalization to 
q-states spins, as in the Potts or Clock model, is straightforward. The problem is to compute the linear response 
R{t, s) on the spin at the site i and at the time t, due to an impulse of external field at an earlier time s and at the 
same site i. Let 



hj{t) = h5,^j9{t - s)9{s + At-t) 



(11) 



be the magnetic field on the z-th site acting during the time interval [s, s-f At], where 9 is the Heavyside step function. 
The response function then is given by [ij, 



R(t, s) = lim 



1 d{a,it)) 



At^o At dhj{s) 



h=0 



where 



dhj{s) 



/i=0 



2^ ^^v{.m,^W I s + At) — 



h=0 



(12) 



(13) 



and [cr] are spin configurations. 

Let us concentrate on the factor containing the conditional probability in the presence of the external field p''([ct'], s-f 
At|[cr"], s). In general, the conditional probability for At sufficiently small is given by 



p{[a],t + At\[aXt) = 5,,,,,n+w{[a'] ^ [a])At + 0{At'), 



(14) 



where we have used the boundary condition p([(t], t| [ct'], t) = Sy^^ y^iy Furthermore, the transition rates must verify 
detailed balance 



w{[a] [a']) exp(-HM/r) = w{W] ^ [a]) cxp(-H[a']/T^), (15) 

where 7i[cr] is the Hamiltonian of the system. 

Introducing the perturbing field as an extra term A7i[(T] = in the Hamiltonian, to linear order in h the most 

general form of the perturbed transition rates w''([cr] — s- [cr']) compatible with the detailed balance condition is 

w\[a] ^ [a']) w\[a] ^ [a']) |l - ^h,{a, ~ a]) + M{[a]., M)| , (16) 

where M([(t], [tr']) is an arbitrary function of order h/T symmetric with respect to the exchange [cr] [cr'], and 
i(;°([cr] [cr']) are unspecified unperturbed transition rates, which satisfy detailed balance. In the following, for 
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simplicity, we shall take M([(t], [cr']) = 0. Implication of this choice, which corresponds to a specification of the 
perturbed transition rates, are discussed in Q. 

Inserting Eqs. ([T4| . (|16p in Eq. (fTS]) . and using the time translational invariance of the conditional probability 
p{[(j], t\[(j'], s + At) = p{[(7], t — At\[a'], s), after some manipulations the response function can be written as 

TR{t,s)^l^^^-lE{t,s) (17) 

where 

Cit,s)^{a,it)a,is)) (18) 

is the autocorrelation function, 

Eit,s)^{a,it)B,is)), (19) 

and 

B. = -J2i<y.~a'J)w"i[a]^[a"]). (20) 



For the dynamic susceptibility one has 



l.„, . „, 1 



Tx{t,s) = -[C{t,t)-C{t,s)]~-J E{t,t')dt', (21) 

It is interesting to observe that Eq. (flTl) is completely analogous to Eqs. ^ and (fTO|) . In fact, it can be easily 
shown that 

'-^-Wt)), (22) 



and that 



^^^^''^ - {B,{t)a,{s)) ^ 0. (23) 



dt 

Subtracting this from Eq. (I17p we finally arrive at Eq. ^ where A{t, s) is given by 



Ait, s) = i [{a,{t)B,{s)) ~ {B,{t)a,is))] . (24) 



Eqs. ([8]) and (I24|) are the main result of this Section. They are identical to Eqs. ^ and (|T0|) for Langevin dynamics, 
since the observable B entering in the asymmetries (jlOp and plays the same role in the two cases. In fact, Eq. (|^^ 
is the analogous of 

'-^-{Bmm) (25) 

obtained from Eq. ([6|) after averaging over the noise. 

In summary, Eq. ([8|) is a relation between the response function and correlation functions of the unperturbed 
kinetics, which generalizes the FDT. Eq. ^ applies to a wide class of systems: Besides being obeyed by soft and 
hard spins, it holds both for COP and NCOP dynamics. Moreover, as it is clear by its derivation, it does not require 
any particular assumption on the Hamiltonian nor on the form of the unperturbed transition rates, and can be easily 
generalized ^1^ to intrinsically non-equilibrium systems where the transition rates do not obey detailed balance. 
Finally, let us briefly discuss (for details see Ref. |8i|) the differences between the results discussed insofar and those 
obtained by Chatelain 12], Ricci-Tersenghi jl^, Diezemann [l^ and Crisanti and Ritort [l3[. Also in these papers, 
in fact, the response function is related to unperturbed correlation functions but, differently from those appairing 
in Eqs. (|17|2ip . these functions must be computed on a system which evolves with an ad hoc kinetic rule, different 
from that of the original unperturbed system, which is introduced with the sole purpose of evaluating the response 
function. It can be shown that this corresponds, in the averaging procedure, to consider only a subset of trajectories 
of the original unperturbed system. Therefore, although the results of Refs. [H, 0, EE Ell important, both for 
computational and analytical calculations, they cannot be regarded as generalizations of the FDT in the sense of 
Eq. ([5]) because the response function is not related to correlation functions of the unperturbed system. 
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IV. FLUCTUATION DISSIPATION RELATION 



In the previous Section we have shown that in the cases considered the integrated response function out of equi- 
hbrium is not only related to the autocorrelation function but also to the correlation E{t, s) by means of Eq. (|21D . 
A very useful tool for the study of slow relaxation phenomena has been introduced by Cugliandolo and Kurchan 
through the off- equilibrium fluctuation dissipation relation (FDR) . This was introduced as a direct relation between 
x{t, s) and C{t, s) as follows: Given that C{t, s) is a monotonously decreasing function of t, for fixed s it is possible 
to invert it and write 



x{t,s)^x{Cit,s),s). (26) 
Then, if for a fixed value of C{t, s) there exists the limit 

hm x{C, s) = S{C) (27) 

s — ^oo 

the function S{C) gives the fluctuation dissipation relation. In the particular case of equilibrium dynamics, FDT is 
recovered and S{C) = [C(t = 0) — C]/T. Originally introduced in the study of the low temperature phase of spin glas s 
mean- field models, the fluctuation dissipation relation has been found in many other instances of slow relaxation [l3|. 



V. STATICS FROM DYNAMICS 



One of the main reasons of interest in the fluctuation dissipation relation is that it may provide a link between 
dynamic and static properties, and in particular with the equilibrium overlap probability function 



^(9) = ^ E exp|-i[i/(M)+i/([a'])]|<5(g([a],[a'])-g) 



(28) 



where Z is the partition function and Q([(7], [cr']) = ^/N'^^aicr'^ is the overlap between two configurations [cr] and 
[a'] . For slowly relaxing systems this is established in general by a theorem by Franz et al. stating that 

1. if S'(C) exists 

2. if limt_>oo x(i, s) — Xeq, Xeq being the equilibrium susceptibility 

then the off-equilibrium fluctuation dissipation relation can be connected to equilibrium properties through 



-T- 



= Piq), (29) 

C=q 



where P{q) is the overlap probability function in the equilibrium state obtained in the limit in which the perturbation 
responsible of x(i, s) is made to vanish. The relation between P{q) and the unperturbed overlap function P{q) must 
be considered carefully. This implies the notion of stocastic stability In a stochastically stable system the 

equilibrium state in the presence of a perturbation, in the limit of a vanishing perturbation, is the same as that of the 
corresponding unperturbed system. Notice that, while stochastic stability is always expected for ergodic systems, this 
property is far from being trivial when more ergodic components are present, as it is easily understood by considering 
the Ising model perturbed by an external magnetic field. If a system is stochastically stable then P{q) — P{q)- A 
milder statement of stochastic stability is that P{q) coincides with P{q) up to the effects of a global symmetry which 
might be removed by the perturbation. For instance, in the Ising case, where the perturbation breaks the up-down 
symmetry, defining 

P{q) = 20{q)P{q) (30) 

the system is stochastically stable in the sense that P{q) — P{q)- In conclusion, if the system is stochastically stable 
Eq. ((29)) holds with P{q) on the right hand side, establishing a connection between the FDR and the equilibrium 
properties of the unperturbed state. On the other hand, if the system is not stochastically stable, P{q) is not related 
neither to P{q) nor to P{q)- As we shall see in Sec. IVI C 21 this is the case of the mean spherical model. 
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With this hnk between statics and dynamics one can translat e llOll to the dynamics the usual classification of 
complex systems based on the kind of replica symmetry breaking [l8|. According to this categorization a first class 
of systems are those whose low temperature phase is characterized by two pure state which are related by a global 
spin inversion. As will be discussed in Sec. lVI Cl these systems without replica symmetry breaking are described by a 
P{q) with a single (5-function centered on the Edwards- Anderson order parameter qea (the squared magnetization, in 
ferromagnetic systems), and their FDR, according to Eq. ([29|) is a broken line with an horizontal part. This situation 
is shown in Fig. [1] upper part (I). A second class of system are those where a transition with a single step of replica 
symmetry breaking occurs, as p-spins with p > 2 in mean field, binary mixtures of soft spheres |19J or Lennard- Jones 
mixtures [20| . In these systems P{q) is made of two (5-functions, one centered in the origin and the other around a 
finite qea- Their FDR is made of two straight lines with finite slopes, as shown in Fig. [T]in the central panel (II). 
Systems as the from zero 
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FIG. 1: Classification of systems according to their P{q) (left) or, equivalently, on S{C) (right), following ref. [K 
represent 5 functions. Rows denoted as I, II, and III describe the three classes of systems discussed in the text. 
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VI. PHASE ORDERING 



Phase ordering [Tl| is usually regarded as the simplest instance of slow relaxation, where concepts like scaling 
and aging, which are the hallmarks of glassy behavior [S^l, can be more easily investigated. However, next to the 
similarities there are also fundamental differences [2l| which require to keep phase ordering well distinct from the out 
of equilibrium behavior in glassy systems, both disordered and non disordered. The main source of the differences 
is the simplicity of the free energy landscape in the case of phase ordering compared to the complexity underlying 
glassy behavior. 

Besides the obvious motivation that the basic, paradigmatic cases need to be thoroughly understood, an additional 
reason for studying phase ordering, among others, is that in some cases the existence of complex slow relaxation 
is identified through the exclusion of coarsening. An example comes from the long standing controversy about the 
nature of the low temperature phase of finite dimensional spin glasses. One argument in favor of replica symmetry 
breaking is that the observed behavior of the response function is incompatible with coarsening ^, llo| . This might 
well be the case; however, for the argument to be sound, the understanding of the out of equilibrium behavior of the 
response function during phase ordering needs to be up to the level that such a delicate issue demands. 
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In this Section we present an overview of the accurate investigation of the response function in phase ordering 
that we have carried out in the last few years. Focusing on the integrated response function ([3]), or zero field cooled 
magnetization (ZFC) in the language of magnetic systems, it will be argued that the response function in phase 
ordering systems is not as trivial as it is believed to be and, after all, it is not the quantity best suited to highlight the 
differences between systems with and without replica symmetry breaking. In fact, as discussed in Sees. IVI CIVI C 21 
there are cases in which phase ordering, and therefore a replica symmetric low temperature state, are compatible with 
a non trivial ZFC. When this happens there is no connection between static and dynamic properties. Phase ordering 
systems offer examples of two distinct mechanism for the lack of this important feature of slow relaxing systems, 
stochastic instability and the vanishing of the scaling exponent of ZFC. 

Let us first briefly recall the main features of a phase ordering process. Consider a system, like a ferromagnet, with 
order parameter (vector or scalar, continuous or discrete) (j){x) and Hamiltonian 7i[(/)(x)] such that below the critical 
temperature Tc the structure of the equilibrium state is simple. For example, in the scalar case, there are two pure 
ordered states connected by inversion symmetry. The form of the Hamiltonian can be taken the simplest compatible 
with such a structure, like Ginzburg-Landau- Wilson for continuous spins or the nearest neighbors Ising Hamiltonian 
for discrete spins. 

Let us generalize the definition ^ to the space and time dependent correlation function 

C{r,t,s) = (0(f,t)0(f' ~ (31) 

where the average is taken over initial condition and thermal noise, and f ~ x — x'. We use the notation C(r — 
0,t,s) = C{t,s), and similarly for the response functions defined below. The linear response function conjugated to 
C(f, t, s) is given by 



(32) 

Sh{x',s) 

where h(x, t) is a space-time dependent external magnetic field and the integrated response function is defined by 

X(r,i,s) = / dsi?(r,t,s). (33) 



A. Dynamics over phase space: equilibration versus falling out of equilibrium 

For a temperature T below Tp, in the thermodynamic limit, the phase space f2 = {[(/)(a;)]} may be regarded as 
the union of three ergodic components [1^ Vl — il+ U il- U fig, where fi± and are the subsets of configurations 
with magnetization hmy^oo y Jy dx(j){x) positive, negative and vanishing, respectively. Denoting by p± [<j){x)\ the 
two broken symmetry pure states, whose typical configurations are schematically represented in Fig. [21 all equilibrium 
states are the convex linear combinations of p±. In particular, the Gibbs state is the symmetric mixture pg[4>{x)\ = 
^ exp(— 7i [(/)(a;)]/r) = ^p^[(f){x)\ + \p^[cf){x)\. The ri± components are the domains of attraction of the pure states 
with p_|_(f2+) = p_(ri_) = 1 and ilo is the border in between them, with zero measure in any of the equilibrium states. 

When ergodicity is broken, quite different behaviors may arise (23| depending on the initial condition pf)[(j){xy\ = 
p{\4){x)\,t = 0). Here, we consider the three cases relevant for what follows, assuming that there are not explicit 
symmetry breaking terms in the equation of motion: 

1. equilibration to a pure state 

if po(f^+) = 1 or /9o(f^-) = 1, in the time evolution configurations are sampled from either one of il± and 
p{[<f>{x)\,t) equilibrates to the time independent pure state p±[(/)(x)] within the finite relaxation time t^q ~ 
where ^ is the equilibrium correlation length and z is the dynamic exponent. The correlation function is the 
same in the two ergodic components and, after equilibration, is time translation invariant 



Cst{r, t) = ((/)(f , i)0(f' , s))± - M2 (34) 

where {(j){x))± = ±M is the spontaneous magnetization. For large distances r ^ ^ and time separations 
t — s ^ teq, the clustering property ((/)(z, i)0(a;', s))± — {4>{x,t))±{(l){x' , s))± is obeyed and the correlations 
decay to zero, as required by ergodicity (see Fig.[S]). 




FIG. 2: Typical configurations of a binary system after equilibration to the pure states p+[(j){x)] or (left and right 

panel) . 



2. equilibration to the Gibbs state 

if /9o(^+) — Poi^-) = 1/2: then configurations are sampled evenly from both disjoint components $7+ and 
The probability density p{[(f>{x)],t) equilibrates now to the Gibbs state paifl'ix)] with the same relaxation time 
teq as in the relaxation to the pure states. Broken ergodicity shows up in the large distance and in the large 
time properties of the correlation function. After equilibration, one has 



CG{r,r) = Cst{f,r) + M^ (35) 
from which follows that correlations do not vanish asymptotically or that the clustering property is not obeyed 



lim CG(f,T) = lim CG(f,T) = M^. (36) 

r — >oo r — *oo 

3. falling out of equilibrium over the border HI . [^ / 

If po{^lo) = 1, for the infinite system p{^lo,t) = 1 also at any finite time after the quench. Namely, the system 
does not equilibrate since in any equilibrium state the measure of fio vanishes. Phase ordering corresponds 
to this case. In fact, the system is initially prepared in equilibrium at very high temperature (for simplicity 
Tj = oo) and at the time t = is suddenly quenched to a final temperature T below Tc- In the initial state 
the probability measure over phase space is uniform po[(j}(x)] — l/|f7|, implying that the initial configuration at 
t = belongs almost certainly to flo, since with a flat measure jfiol is overwhelmingly larger than \ fl±\. 

The morphology of typical configurations visited as the system moves over fig is a patchwork of domains of the 
two competing equilibrium phases, which coarsen as the time goes on, as schematically shown in Fig. [3l The 
typical size of domains grows with the power law L{t) ^ t^^^ , where z = 2 (independent of dimensionality) for 
dynamics with non conserved order parameter ^TJ], as it will be considered here. The sampling of configurations 
of this type is responsible of the peculiar features of phase ordering. At a given time s there remains defined a 
length L{s) such that for space separations r ^ L(s) or for time separations t — s <^ s intra-domains properties 
are probed. Then, everything goes as in the case 2 of the equilibration to the Gibbs state, ergodicity looks broken 
and the correlation function obeys Eq. p5p . Conversely, for r ^ L{s) or t/s S> 1, inter-domains properties are 
probed, ergodicity is restored (as it should be, since evolution takes place within the single ergodic component 
fio) and eventually the correlation function decays to zero. However, the peculiarity is that if the limit s oo 
is taken before r — > oo, in the space sector ergodicity remains broken giving rise, for instance, to the growth of 
the Bragg peak in the equal time structure factor. 

According to this picture, the correlation function can be written as the sum of two contributions 



C(f, t, s) = Cst{r, t) + Cag{r, t, s) 



(37) 






FIG. 3: Configurations of a coarsening system at different times ti < t2 < < t4. 



where the first one is the stationary contribution of Eq. (|34p describing equihbrium fluctuations in the pure 
states and the second one contains all the out of equilibrium information. The latter one is the correlation 
function of interest in the theory of phase ordering where, in order to isolate it, zero temperature quenches are 
usually considered as a device to eliminate the stationary component. It is now well established that C'ag{r, t, s) 
obeys scaling in the form ^26i | 



Cag{r,t,s)^C{r/L{s),t/s) 



(38) 



with C(x, y) — M'^ for a; < 1 and y ^ 1, while 



C{r/L{s),t/s)^{t/s)-^/^-h{r/Lis)) 
for large time separation [Tlj . where A is the Fisher-Huse exponent. 



(39) 



B. Zero field cooled magnetization 

Let us next consider what happens when a time independent external field h{x, s) is switched on at the time s. To 
linear order the expectation value of the order parameter at the later time t is given by 



(0(f , t))h = (0(f , t))o + J dx' J dsR{x - x', t, s)h{x', s) (40) 

If h[x^ s) — h{x)9{t — s) is a random field switched on and kept constant from s onwards, with expectations 

/i(f) = (41) 



h{x)h{x') = h^5{x-x') 



then one has 



X{x-y,t,s) = h {(j){x,t))hh{y). 



(42) 
(43) 
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Namely, ZFC is the correlation at the time t of the order parameter with the random external field. 

Going to the three processes considered above, and restricting attention from now on, for simplicity, to the case of 
coincident points (r — 0) 

1. after equilibration in the pure state has occurred and the stationary regime has been entered, the order parameter 
correlates with the external field via the equilibrium thermal fluctuations, FDT is obeyed 

Xst{T) = ^[CstiT = 0)-Cst{T)] (44) 

and since Cst{T) decays to zero for t > t^q, over the same time scale Xst{T) saturates to 



lim Xstir) = Xeq = Tf^Ceq (45) 

which is the susceptibility computed in the final equilibrium state p±[(j){xy\ (see Fig. [6]). 

2. As far as ZFC is concerned, there is no difference between the relaxation to the mixed Gibbs state and the 
relaxation to a pure state. Hence, FDT is satisfied and can be written both in terms of Cst or Cg since, as 
Eq. ((35|) shows, they differ by a constant. 

3. In the phase ordering process the system stays out of equilibrium, so it useful to write ZFC as the sum of two 
contributions ^] 

X{t,s)^ Xstir) +Xag{t,s) (46) 

where Xstir) satisfies Eq. ((44l) and Xag{t, s) represents the additional out of equilibrium response. In connection 
with this latter contribution there are two basic questions 

i) how does it behave with time 

ii) what is the relation between Xag and Cag, if any. 

1. Scaling hypothesis 

Since ZFC measures the growth of correlation between the order parameter and the external field, the first question 
raised above addresses the problem of an out of equilibrium mechanism for this correlation, in addition to the thermal 
fiuctuations accounting for Xst- The starting point for the answer is the assumption of a scaling form 

Xag{t,s)^ S-''-Xag{t/s) (47) 

which is the counterpart of Eq. ([38)1 for the correlation function. 

The next step is to make statements on the exponent a^^^ and on the scaling function Xag{x) . There exists in the 
literature an estimate of based on simple reasoning. What makes phase ordering different from relaxation in the 
pure or in the Gibbs state is the existence of defects. The simplest assumption is that Xag{t, s) is proportional to the 
density of defects dllSlli. This imphes 

= 6 (48) 
where the exponent 6 regulates the time dependence of the density of defects p{t) ^ L{t)^^ ^ t^^ , namely 



S — n/z 

with n = 1 for scalar and n — 2 for vector order parameter [ll| 



(49) 
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According to this argument should be independent of dimensionahty. This conclusion is not corroborated by 
the available exact, approximate and numerical results. On the basis of exact analytical solutions for the d = 1 Ising 
model [2^ [lOI and for the large N model [sH , approximate analytical results based on the Gaussian auxiliary field 
(GAF) approximation [H, [s^ and numerical results from simulations [H, [s^, H^, [H, [13] with d = 2,3,4, we have 
argued that 



'^x ~ { S with log corrections for d ~ du (^'^) 
[ (5 for d> du 

where and du > dL do depend on the system in the following way 

• (ii is the dimensionality where = 0. In the Ising model = 1, while in the large N model d^ = 2. The 
speculation is that in general dL = 1 for systems with discrete symmetry and = 2 for systems with continuous 
symmetry, therefore suggesting that coincides with the lower critical dimensionality of equilibrium critical 
phenomena, although the reasons for this identification are far from clear. 

• d[/ is a value of the dimensionality specific of ZFC and separating d < du, where depends on d, from d > du 
where is independent of dimensionality and Eq. (|48p holds. The existence of du is due [sHj to a mechanism, 
i.e. the existence of a dangerous irrelevant variable, quite similar (including logarithmic corrections) to the one 
leading to the breaking of hyperscaling above the upper critical dimensionality in static critical phenomena. 
However, du cannot be identified with the upper critical dimensionality since we have found, so far, d;/ = 3 in 
the Ising model and d^/ = 4 in the large N model. In the scalar case it may be argued [36, 38] that du coincides 
with the dimensionality djj = 3 such that interfaces do roughen for d < d^, and do not for d > du- This will be 
discussed in Sec. IVIB 21 A general criterion for establishing the value of du, however, is not yet known. 

The validity of Eq. (|47)) with O;^ given by Eq. (|50p has been checked, in addition to the cases where analytical results 
are available, with very good accuracy in the simulations of the Ising and clock model and of the time dependent 
Ginzburg-Landau equation [H, [13, [HI, [s^, [s^- The values of 5, dL and du obtained for the different systems are 
collected in Table |T] and the behavior of as dimensionality is varied is displayed in Figd) 





Ising 


GAF 


N = 00 


6 


1/2 


1/2 


1 


dh 


1 


1 


2 


dx 


3 


2 


4 



TABLE I; Parameters entering Ea.(|50p in various models. 



2. Roughening of interfaces 

Apart from the few exact solutions mentioned above there is not a general derivation of Eq. ([50)) which, at this stage 
remains a phenomenological formula. For the case of a scalar order parameter, an argument has been proposed [36l[38j 
explaining the dependence of on d in terms of the roughening properties of the interfaces. It is based on two simple 
physical ingredients: a) the aging response is given by the density of defects p{t) times the response of a single 
defect [13 Xag{t,s) — p{t)Xag(.^T ^) ^) each defect responds to the perturbation by optimizing its position with 
respect to the external field in a quasi-equilibrium way. In d = 1 this occurs via a displacement of the defect [ssj . 
In higher dimensions, since defects are spatially extended, the response is produced by a deformation of the defect 
shape. 

We develop the argument for a 2-d system, the extension to arbitrary d being straightforward. A defect 
is a sharp interface separating two domains of opposite magnetization. In order to analyze x^g(t, s) we con- 
sider configurations with a single defect as depicted in Fig. [5l The corresponding integrated response function 
reads @ X^giti^) = l/{h'^C'^~'^) f dxdy ((/)(x, y))/i(x, y), where 0(x, y) is the order parameter field which satu- 
rates to ±1 in the bulk of domains, and x, y are space coordinates. /i(x, y) is the external random field with 
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□ 



NCOP N = 1 (expected) 

COP N = 1 (expected) 

NCOP N > 1 (expected) 

COP N > 1 (expected) 

NCOP N = 1 Numerical 

COP N = 1 Numerical 



0,6 <^ NCOP N > 1 Numerical 



N > 1 Numerical 




FIG. 4: Exponent in different coarsening systems with scalar and vector order parameter, non-conserved and conserved 
order parameter The continuous lines represent Eq. (I50|) . while the dots are the values from simulations [ssl. [3^ . 




FIG. 5: Configurations with a single interface at time s (dashed line) and at time t (continuous line). 



expectations (|41|42p . and C is the linear system size. The overbar and angular brackets denote averages over 
the random field and thermal histories, respectively. With an interface of shape Zs{y) at time s (Fig. [5]), we 
can write Xagit^^) = S/i P;i({z(y)}, i), where Ph{{z{y)},t) is the probability that an interface 

profile {^(y)} occurs at time t and Eh = — dy J^^^^^ dxh{x, y)sign[z{y) — Zs(y)] is the magnetic energy. We 
now introduce assumption b) making the ansatz for the correction to the unperturbed probability Po{{z{y)},t) 
in the form of a Boltzmann factor P,,({z(y)}, i) = PQ{{z{y)},t) exp{-Eh/T) ~ Poi{z{y)},t)[l - Eh/T]. Then 

Xigitys) = -E;i(l — -E/i/r)Po({^(y)}, 0- Taking into account that the term Hnear in Eh van- 

ishes by symmetry and neglecting Zgij) with respect to z{y) for i ^ s, we eventually find TXagit, s) = 
J^^j dy|z(y)|Po({z(y)}, i). This defines a length which scales as the roughness of the interface given by 
W{t) = [C^-'^J^^y J (iyz(y)2Po({2(y)},t)]i/2. The behavior of W{t) in the coarsening process can be inferred from an 

argument due to Villain |39[. In the case d < 3, when interfaces are rough [40|, for NCOP one has W{t) ~ t(3-'^)/4, 
while for COP W{t) ~ t(3-<i)/6^ with logarithmic corrections in both cases for d = 3. For d > 3 interfaces are flat 
and W{t) ~ const. Finally, multiplying xlg by p{t) ^ L{t)^^ Eq. (|5D|) is recovered [4l| and du is identified with the 
roughening dimensionality du = 3. 



C. Statics from dynamics 



We may now check if, and how, the connection between statics and dynamics discussed in Sec. |V] is realized in 
phase ordering systems. In the following we shall consider d > d^. 

In order to search for S{C) in the case of phase ordering, let us set r = in Eq. ([55)1 and let us eliminate t/ s 
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between Xag and Cag obtaining 

Xagit,s)^S-'^-XagiCag). (51) 

Then, from Eqs. (|46l44l5ip one can write the general relation 

Xit, 5) = ^ [Cstir = 0) - C^tir)] + S-'^-XagiCag). (52) 

Using the identity [Cst{T = 0) — Cst{T)] — [Cst{T = 0) + Af^ — Cst{T) — Af'^] and considering that, as shown 
schematically in Fig.[6l in the time interval where Cst{T) ^ 0, i.e. for short times, one can replace Cag{t/s) with A'P 
or equivalently Cst{T) + M"^ = C{t, s), the above equation can be rewritten as 

X{t,s)^Xst{C)+S~''-XagiCag) (53) 

where the function XstiC) is defined by 

_ / [C{t, t) - C{t, s)] for NP<C< C{t, t) 

j [C7(t,t)-Af2] for C<M\ ^^^^ 

Therefore, from Eq. (|53p we have that for phase ordering systems the fluctuation dissipation relation exists if a^^ > 
(i.e. for d > d^) and it is given by 

S{C)^Xst{C). (55) 
Computing the derivative in the left hand side of Eq. and using Eqs. ([55)1 and for d > we find 



dC^ 



= 5{q - A.f2). (56) 

C=q 



Coming to statics, in replica symmetric low temperature states, as for instance in ferromagnetic systems, the overlap 
function is always trivial and, as anticipated in Sec. |Vl one has 

m = \ W - M^) + 5{q + Af 2)] , (57) 

with 

P{q) = P{q) = S{q - AP). (58) 

as shown in Fig. [1](I) (we recall that qEA = in this case). From Eqs. (|58l56p . therefore, Eq. (P^)) is satisfied, and 
the connection between statics and dynamics holds. 

For flpj, = a little more care is needed. Equation (|55|) yields xi^i^) = Xst{C) + Xag[Cag)- Recalhng that — Q 
occurs at d = c?l, which coincides with the lower critical dimensionality, in order to have a phase ordering process a 
quench to T = is required. This, in turn, implies Cst{t, s) = and Cag{t, s) — C{t, s). Therefore, using Eq. ([54]) we 
have 



SiCj^Xlg+XagiC) (59) 

where xlq = linir— >o [C'(O) ~ Af^J/T is the T = equilibrium susceptibility, which vanishes for hard spins while is 
different from zero for soft spins. Therefore the FDR exists also in this case. However, while for > Xag eventually 
disappears and Eq. (|45p holds, this is no longer true for d = dL- Here = Q and, consequently, as can be seen from 
Eq. (|59)) Xag gives a contribution to the response which persists also in the asymptotic time region. Then Eq. (|45|) . 
and hence condition (2) above Eq. (^51 are not fulfilled. Being one of the hypothesis leading to Eq. violated, 
the connection between statics and dynamics could not hold. Actually, in all the model explicitly considered in the 
literature [s^, [H, [s^ it turns out that at d = S{C) is a non-trivial dynamical function unrelated to P{q)- For the 
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FIG. 6; Schematic plot of the behavior of the two time functions in the coarsening stage as f — s is varied keeping s fixed 
(upper figure), and of the resuhing fiuctuation dissipation plot (lower figure). The stationary parts Cst, Xst, as discussed in 
Sees. IVI AlVI Bl saturate to their final value on times r ~ tsq; the aging parts Cag, Xag, according to the scaling forms (|38I47|) 
remain constant up to times t — s ~ s. The magnitude of Xag is proportional to and decreases as s is increased (shown 

by an arrow in Figure). 



sake of definiteness, let us discuss the case of the Ising model with d = 1 [29|, |30|. In order to make compatible the 
two requirements of having an ordered equilibrium state and a well defined linear response function, instead of taking 
the T ^ limit it is necessary to take the limit of an infinite ferromagnetic coupling [2§|. Then, P{q) and P(q) are 
given by Eqs. ([57)1 and (|55|) with = 1 at all temperatures. On the other hand, for any T one also have Q (see 
Fig. ED ' 



V2 

Txag{C) = — arctan 



\/2cot (I 



C 



(60) 



This gives 



C=q 



TT cos(7rg/2) sin(7rg/2) 
[2 - sin(7rg/2)]2 ' 



(61) 



Hence, it is clear that Eq. (|29p is not verified. The reason is that the second of the above conditions required for 
establishing the connection is not satisfied. In fact, from Eqs. (j53p and (|60p . keeping in mind that the limits t ^ oo 
and C are equivalent, we have 



t- 



lim Txit, s) = 1/V2 (62) 



which is responsible of the violation of condition (2) above Eq. (|29p . since in this case Xeq = 0. Interestingly, a 
similar behavior is observed p^j also for the Ising model on graphs with Tc = 0, which, in a sense, can be regarded 
as being at di,. 



1. Role of quenched disorder at d = dL 



The behavior of the exponent a^, and its vanishing at = 1 can be qualitatively interpreted in terms of the behavior 
of the response XagC^i •s) associated to a single interface. In d = 1 it can be shown exactly [sl] that Xag(i, s) oc (i— s)^/^. 
Therefore, when computing the total response through Xag{t, s) ~ p{t)x''{t, s) the loss of interfaces described by p{t) 
is exactly balanced by the increase of x*(^j s)j which leads to a finite limt^oo Xag{t, s). This in turn is responsible for 
the breakdown of the condition (2) above Eq. (|^ . For d > dL, instead, the growth of x'*(ij s) is not sufficient ^] to 
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FIG. 7: Fluctuation dissipation plot for the Id-Ising model. 



balance the decrease of p(t). This happens because, while in c? = 1 interfaces are Brownian walkers, free to move in 
order to maximize the response, for d > 1 this issue is contrasted by surface tension, restoring the validity of condition 
(2). In this Section how a similar effect, namely a reduction of the response of interfaces, can be produced, also in 
d = 1, by the presence of quenched disorder. 

Let us consider the case of the d = 1 random field Ising model (RFIM). In the presence of a quenched random 
field, domain walls perform random walks in a random potential of the Sinai type and the average domain size L(t) 
behaves as the root mean square displacement of the random walker (43j. The typical potential barrier encountered 
by a walker after traveling a distance I is order of ^/l<Jh where ah is the variance of the random field. Hence, there 
exists a characteristic length Lg ~ T"^ /au representing the distance over which potential barriers are of the order of 
magnitude of thermal energy. For displacements much less than Lg diffusion takes place in a flat landscape like in 
the pure system and L{t) ^ t^l"^ . For displacements much greater than Lg, instead, one finds the Sinaijl^] diffusion 
law L(i) ~ (Ini)^. The response function obeysfi^ the scaling relation (Fig. [S]) 



rXag(i,S,^g 



Lisy 



(63) 



where z = L{s)/Lg. For z = the form of the response function for the pure system is recovered. With z > there 
is a crossover. The pure case behavior holds for L{t) — L{s) <C Lg, while for L{t) — L{s) > Lg the response levels off 
and then decreases. This is clearly displayed also in the plot (Fig. [5]) against the autocorrelation function. Looking at 

the effective response of a single interface x*(<, s) one finds x''(i, s, Lg) = L{s)x'^ (x(s7' •^j "^i^h the scaling function 
displaying the behavior 



X"(x,z) 




0) , for a; - 1 < 1/z 
, for cc — 1 >• 1/z. 



(64) 



From this follows X^i't^^) ^ the preasymptotic regime and x'*(i, s) ~ Pi^^^i^) i^i the asymptotic regime, 

which account for the crossover of the response function in Figs. l8|9l in terms of the balance between the rate of growth 
of the single interface response and the rate of loss of interfaces. Hence, for z > eventually xi^j vanishes and 
in the limit z — s- oo one expects Xag(a^, ^ ~ oo) = 0. Therefore, for any finite quenched random field the validity of 
Eq. (gll) is restored. 
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2=0 




FIG. 8: The ZFC Xaa{t, s) is plotted versus the the ratio [L{t) - L{s)]/L{s), for Lg = 400, 1600, and z = 0.1, 9. The sohd line 
is the exact result for z — > 0. 
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FIG. 9: The ZFC Xagi^, s) is plotted versus the autocorrelation function C{t, s), for Lg — 25, 100, and z = 1,4, 9, 16. The solid 
line is the exact result for z 0. 



2. Failure by stochastic instability 

An interesting example [46l | , where statics cannot be reconstructed from dynamics because the third requirement of 
stochastic stability is not satisfied, comes from the spherical model. More precisely, one must consider in parallel the 



original version of the spherical model (SM) of Berlin and Kac [47[ and the mean spherical model (MSM) introduced 
by Lewis and Wannier (48|, with the spherical constraint treated in the mean. These two models are equivalent above 
but not below Tc (49j . The low temperature states are quite different, with a bimodal order parameter probability 
distribution in the SM case and a Gaussian distribution centered in the origin in the MSM case. The corresponding 
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overlap functions are also very different [46|. Considering, for simplicity, T = one has 




where Kq is a Bessel function of imaginary argument (FigllOp. 




-10 1 q 



FIG. 10: Overlap distribution for mean spherical model with = 1. The arrows represent the S functions of the overlap 
distribution for the spherical model. 



However, after switching of an external field, one finds for both models P{q) = S{q—AP). This means that stochastic 
stability holds for SM but not for MSM. 

On the other hand, the relaxation properties are the same in the two models, both above and below Tc if the 
thermodynamic limit is taken before the t ^ oo limit (46j . Then, the linear response function is the same for both 
models and obeys Eq. with given by Eq. where S, and du are the same as for the large N model 

(Table I). Hence, we have that although Eq. is satisfied for both models, nonetheless statics and dynamics are 
connected only in the SM case, where P(q) = P(q)- Instead, this is not possible in the MSM case where P{q) ^ P{q)- 

VII. CONCLUSIONS 

In this Article we have reviewed some recent progresses in the field of non-equilibrium linear response theory. 
A first accomplishment is the derivation of a generalization of the FDT for Markov processes which allows the 
computation of the response function in terms of correlation functions of the unperturbed system. This represents a 
great simplification particularly in numerical calculations, which are usually computationally very demanding: The 
generalization of the FDT allows a sensible speed and precise numerical determination of the response function can 
be achieved. This quantity has been deeply investigated particularly in the field of slowly relaxing systems, because 
its relation with the autocorrelation function represents a bridge between statics and dynamics. 

Phase-ordering systems can be regarded as the simplest instance of aging systems, where the behavior of the response 
function can be more easily investigated. In this context, a partial understanding has been achieved by matching 
the results of numerical simulations with the outcomes of solvable models and approximate theories, showing that 
the scaling properties of the response function are non-trivial. In particular, Eq. (|T7)l is obeyed with depending 
on dimensionality through the phenomenological formula (j50|) . which is found to be consistent with all the cases 
considered in the literature and, for the scalar case is supported by an argument based on the roughness properties of 
the interfaces. The dependence of on dimensionality is such that it vanishes at the lower critical dimension. This 
implies an asymptotic finite contribution of the aging part of the response function which invalidates the connection 
between statics and dynamics. Phase ordering therefore offers examples where a replica symmetric low temperature 
state is compatible with a non trivial FDR which, therefore, cannot be used to infer the properties of the equilibrium 
state. 

This whole phenomenology is not adequately captured by the existing approaches to phase-ordering. Theories based 
on the GAF method, originally introduced by Otha, Jasnow and Kawasaki [s^l, provide the phenomenological for- 
mula (jSOp but with a wrong value djj = 2 [3^ . [33| . This discrepancy is not removed using a perturbative expansion |5lj 
developed to improve over the GAF approximation. Next to these theories, it is of much interest the approach by 
Henkel et al. |52l |. based on the conjecture that the response function transforms covariantly under the group of local 
scale transformations. This ansatz, however, fixes the form of the scaling function in Eq. (j47p but not the exponent 
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which remains insofar an undetermined quantity. A first principle theory for the complete description of the behavior 
of the linear response function in phase-ordering systems may represent a pre-requisite for understanding the behavior 
of more complex systems, like glasses and spin glasses. However, despite some progresses of a specific character, such 
a theory is presently still lacking. 
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